Muonium as a shallow center in GaN 
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A paramagnetic muonium (Mu) state with an extremely small hyperfine parameter was observed 
for the first time in single-crystalline GaN below 25 K. It has a highly anisotropic hyperfine structure 
with axial symmetry along the (0001) direction, suggesting that it is located either at a nitrogen- 
antibonding or a bond-centered site oriented parallel to the c-axis. Its small ionization energy (< 14 
meV) and small hyperfine parameter (~ 10"** times the vacuum value) indicate that muonium in 
one of its possible sites produces a shallow state, raising the possibility that the analogous hydrogen 
center could be a source of n-type conductivity in as-grown GaN. 

PACS numbers: 61.72.Vv, 71.55.Gs, 76.75.+i 
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Since the discovery of methods to produce sufficient 
p-type conductivity by Mg-doping gallium ni- 

tride and related compound semiconductors are being ag- 
gressively developed for electronic and optelectronic de- 
vices such as blue/green lasers and light-emitting diodes. 
Unique features such as a wide and direct band gap and 
high breakdown field make the nitrides ideal for such 
applications. However, as-grown undoped GaN epitax- 
ial thin films, as well as bulk single crystals, commonly 
exhibit n-type conductivity with concentrations ranging 
from 10^® to 10^^ cm^'^. Extensive experimental and the- 
oretical studies 3] have been undertaken to understand 
the origin of this n-type conductivity. 

For many years, nitrogen vacancies, which are com- 
monly observed in as-grown GaN, were thought to be a 
major source of the n-type conductivity Recent 
theoretical work challenges this view by showing that 
the nitrogen vacancies have a high formation energy in 
n-type GaN; hence their concentration is predicted to 
be too low to affect the electrical conductivity. More- 
over, Hall effect measurements on a GaN sample irra- 
diated by high-energy electrons showed that the donor 
level of the nitrogen vacancies (64(10) meV) was much 
deeper than that of the residual donors (18 meV) 0- 
Contamination by oxygen or silicon is difficult to avoid 
during crystal growth and both of these impurities have 
been proposed as the origin of unintentional n-type be- 
haviour. Silicon is often used as an intentional donor 
dopant and can provide free electron concentrations of 
up to lO^^cm"^. Very recent magneto-optical studies 
of hydride-vapor-phase epitaxial GaN show that a can- 
didate other than O or Si impurities can also act as a 
shallow donor 

By contrast, magnesium is currently the only acceptor 
dopant reliably used to obtain p-type GaN. Even in this 



case, high hole concentrations were difficult to achieve 
until it was found that hydrogen reacts with and passi- 
vates the Mg acceptors 9] , similar to its effect on accep- 
tors in Si. A post-growth anneal at high temperature is 
required to remove H from the vicinity of the Mg dopants 
and thereby activate p-type electrical properties. 

Motivated by such a crucial role of hydrogen in GaN, 
extensive fiSK (muon spin resonance) studies ^3 have 
been performed to clarify the physical and electronic 
structure of isolated H centers via their muonium ana- 
log. The results established the existence of two charged 
states (Mu+ and Mu") [ll[l3, and level-crossing reso- 
nance spectra revealed that Mu+ and Mu~ reside at sites 
anti-bonding to N and Ga, respectively [isj . 

In addition, neutral paramagnetic muonium states, i.e. 
Mu*^, are readily observed in a wide variety of semicon- 
ductors. While the dynamical aspects (e.g., diffusion) 
may be considerably different between Mu and H due to 
the light mass of Mu (m^ ~ ^Wp), the local electronic 
structure of Mu is virtually equivalent to that of H after a 
small correction due to the difference in the reduced mass 
4%). Recently, a novel Mu state having an extremely 
small hyperfine parameter (10"** x A^) was reported in 
H-VI cornpound semiconductors including CdS 0| and 
ZnO [iEEIH, implying that Mu (and hence H) can 
act as a donor in these materials. No Mu" hyperfine 
spin-precession spectra have previously been reported for 
GaN, although hyperfine decoupling measurements sug- 
gest a short-lived atomic-like neutral 13J. 

In this Letter we report on the first observation of a 
paramagnetic muonium spectrum in GaN. The observed 
Mu° state has an extremely small and highly anisotropic 
hyperfine parameter. The location within the band gap 
for the [0/-f] energy level associated with this Mu state is 
estimated from the measured activation energy for ther- 
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FIG. 1: The /iSR time spectrum in GaN at 5.0 K under an 
external field of 3.0 T applied parallel to the [0001] axis, dis- 
played in a rotating-reference-frame frequency of 405 MHz. A 
beat pattern due to multiple frequencies is clearly seen. 

mal ionization. These results imply that an isolated hy- 
drogen impurity would behave as a shallow donor if it 
were located at the same crystallographic site. 

We performed /iSR measurements on a GaN single- 
crystal with the hexagonal (2H) wurtzite structure. GaN 
has tetrahedral coordination with bonds parallel to the 
symmetry axis slightly elongated and off- axis bond di- 
rections ordered to give hexagonal symmetry around 
[0001]. The /uSR experiment was conducted at the TRI- 
UMF M15 bcamlinc with the Belle spectrometer. Muons 
from a surface beam (100% spin-polarized with an en- 
ergy of 4 MeV) with their polarization transverse to 
the applied magnetic field were implanted into GaN 
(3 X 3 X 0.5 mm, [0001] orientation, n type with a concen- 
tration of 10^^ cm~^). The subsequent time evolution of 
the muon spin precession yields information on the muo- 
nium hyperfine parameters and stability of the center. 
High fields (up to 5 T) were used to quench fast spin re- 
laxation due to interactions with host nuclear spins. The 
GaN crystal was mounted in a He gas-flow cryostat such 
that the [0001] and [1120] crystallographic axes could be 
oriented at specific angles to the applied field in order to 
examine the symmetry of the hyperfine interaction. De- 
pendences on temperature and magnetic field strength 
were obtained with the field parallel to the [0001] axis. 

Above 25 K, only a single (diamagnetic) precession 
signal is observed at the muon Larmor frequency (gy- 
romagnetic ratio 7^ = 27r x 135.53 MHz/T). Relaxation 
of this signal is well described by Gaussian damping with 
a nearly temperature-independent rate constant of ~ 0.2 
;t<s~^. This damping rate is satisfactorily explained by 
the dipole-dipole interaction of muons with ^'•^■''^Ga and 
^^N nuclei. The muon spin rotation signal changes dras- 
tically below 25 K. A typical /iSR spectrum is shown in 
Fig. 1. Fig. 2 shows the angular and temperature de- 
pendence of the frequency spectra obtained by Fourier 
transform, in which one pair of satellite lines is clearly 
seen with their positions situated symmetrically around 
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FIG. 2: Frequency spectra obtained for GaN at (a) 2.5 K and 
(b) 22.5 K with 5=3.0 T parallel to [0001] axis, and with 
5=1. 5T at 2.5 K, where [0001] is parallel to B (c) or tilted 
by 35° from B (d), where from 0.01 /xs to 9.60 /xs is adopted 

as a time range of FFT. The number of y-axis in (a) and the 
other is not directly comparable, because of the difference in 
the data statistics. 

the central line, which corresponds to the precession of 
diamagnetic muons. The splitting of these satellites re- 
mained unchanged when the applied field was increased 
from 1.5 T up to 5 T (see Fig. 2a and 2c), a result that is 
important in identifying the spectra as due to the hyper- 
fine interaction of a Mu° center. The splitting decreases 
when the [0001] axis is tilted with respect to B as in 
Fig. 2d. Moreover, an equivalent frequency spectrum was 
observed when the [1120] axis was rotated by 90° around 
[0001], which was oriented at 35° to the applied mag- 
netic field. These observations demonstrate the presence 
of a paramagnetic muonium state in GaN. The result- 
ing hyperfine interaction is extremely small, about 10"** 
times the vacuum value for a Mu atom, and is axially 
symmetric with respect to [0001]. 

With the [0001] axial symmetry established qualita- 
tively, a more formal analysis can be undertaken to ex- 
tract the hyperfine constants. The amplitude differences 
for the two hyperfine lines allow the specific muonium 
transition to be assigned to each satellite and estab- 
lishes the sign of the hyperfine constants. Spin preces- 
sion signals from the paramagnetic Mu state have two 
components in the high- field limit {B ^ 2TrA/je, where 
7e = 27r X 28.024 GHz/T is the electron gyromagnetic 
ratio): 

,y,2{0) ^1^0- \^y{e), (1) 

v^i{e) + (2) 

Av{e) = A{e)=A^^cos'^e + A_Lsm'^0, (3) 

where 2^1/0 = Ji^B, 6 is the angle between B and the 

[0001] symmetry axis, and j4[j and A±_ are the hyperfine 
parameters parallel and perpendicular to [0001], respec- 
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FIG. 3: Magnetic field dependence of tlie ratio of the lower- 
frequency satellite amplitudes to the sum of both satellites 
at 2.5 K (symbols). Solid curve is a fitting result using the 
Maxwell-Boltzmann distribution with the Mu temperature as 
a free parameter. 

tively. The frequencies ^12 and 1^34 correspond to the 
transitions of spin states between |se,s^) = 1 + 5,+^) 
and I + ^,—5), and between | — ^,+5) and | — 5,-5), 
respectively. 

Because the population of the Se — —■^ state is always 
larger in an applied field, the observations in Fig. 2 im- 
ply that the low-frequency line corresponds to Se = +5, 
thus to 1^12, and that A{0) is positive for the displayed 
orientations. Results from direct fits to the time-domain 
spectra (see Fig. 1), where 0.01 /iS to 9.60 ^s in time 
range, provide a better amplitude measurement for com- 
parison with theoretical expectations. Analysis assum- 
ing three frequency components {vq, and 1^34, we 
assumed the depolarization rates of the two Mu lines are 
equal ) is quite satisfactory, with a reasonably small re- 
duced (— 1-55) for the data in Fig. 1. Fig. 3 shows 
the field dependence for the ratio of the lower-frequency 
satellite amplitudes to the sum of both satellites at 2.5 K 
and 9 = 0, together with a fitting result of electron spin 
occupation probabilities for an isolated muonium center 
unsing the Maxwell-Boltzmann distribution, where 3.4 K 
was obtained as the Mu temperature. Correspondence 
between the sample and obtained Mu temperature is also 
satisfactory, confirming assignment of a positive hy- 
perfine constant. 

From the spectrum in Fig. 2a, with B applied along 
[0001], the hyperfine constant is deduced to be 

A{0°) = A|| = +337(10) kHz. (4) 

Combining this result with the data for a tilted sample 
orientation with respect to B (Fig. 2, where 9 = 35.0°, 
and Av = 146(3)kHz) the remaining hyperfine constant 
is found to be 

= -243(30) kHz. (5) 

As an experimental check, these parameters reproduce 
the observed splitting for an additional orientation with 
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FIG. 4: The fractional yield of Mu (closed circles), and a 
diamagnetic state (open circles) versus temperature in GaN. 
Solid curves are fitting results (see text). 

9 = 15.0°. Therefore, we conclude that the hyperfine 
parameters parallel and perpendicular to the [0001] sym- 
metry axis have opposite signs. The static dielectric con- 
stants in GaN are 10.4 parallel and 9.5 perpendicular to 
[0001] 3]. Since the degree of anisotropy in the hyperfine 
tensor for the observed paramagnetic muonium state is 
much larger than that implied by the dielectric constant, 
we conclude that the observed anisotropy is due to lo- 
cal site symmetry. Of the likely locations for Mu in the 
GaN wurtzite structure [lallSlliflli these results strongly 
suggest that the paramagnetic Mu center resides at either 
the ABjvii or the BCy site. 

The temperature dependence of the amplitudes of the 
Mu° and diamagnetic signals are plotted in Fig. 4. The 
total amplitude is almost independent of temperature, 
suggesting that Mu is ionized to a diamagnetic state 
above 25 K. It seems unlikely that the energy levels are 
just above the valence band since, if that were correct, 
Mu should remain neutral due to the bulk n-type conduc- 
tivity of the present specimen which puts the Fermi level 
much higher than mid-gap. These results indicate that 
the observed paramagnetic Mu center acts as a shallow 
donor. Thus, since Mu centers simulate the electronic 
structure of H in GaN, our result provides convincing 
evidence that hydrogen centers at this particular site are 
shallow donors, and thus would contribute to n-type con- 
ductivity in GaN. 

The activation energy Ea for Mu ionization was found 
to be 4.6(6) meV by fitting the data in Fig. 4 over a region 
from 10 K to 25 K with a function a + (3exp{—Ea/kBT). 
Interpretation of this energy depends on the precise pro- 
cess forming the neutral charge-state of the Mu donor and 
the extent to which equilibrium arguments apply to muo- 
nium state occupations. At one end of the range it repre- 
sents a direct transition from the defect level to the bot- 
tom of the conduction band; while at the other extreme. 
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in equilibrium it is for a transition from the defect level to 
the Fermi energy. Assuming minimum compensation, a 
reasonable estimate of Ep at non-zero temperatures be- 
low the ionization of residual donors is midway between 
the 18 meV level for the dominant donors and the bottom 
of the conduction band. Thus our ~5 meV activation en- 
ergy translates to a Mu[0/-|-] level somewhere between 5 
and 14 meV below the conduction band edge. Consider- 
ing the difficulty in interpreting the Mu ionization energy 
and ambiguity in determining the level of the shallower 
low-concentration donor in Ref. the hydrogen equiv- 
alent of the observed Mu° center might be a candidate 
for this additional donor. 

It is interesting that for Mu/H in GaN, no such shal- 
low states are predicted in the latest theoretical studies 
[20I EH . Our observation may not be incompatible with 
those predictions if the muon location that yields the 
shallow center were metastable rather than the lowest en- 
ergy site. Of the two sites that locally satisfy the hyper- 
fine symmetry, BC|| is one of the lower-energy positions 
for both neutral and positive change states. Data for hy- 
drogen at high temperatures indicate that both the 
ABjv± and BC|| sites are occupied during H+ diffusion. 
The other likely site, ABtvh, is at a significantly higher 
energy, but yields a (meta)stable location for Mu+ based 
on level-crossing results (UJ. The same 'cage' region of 
the GaN structure that contains ABjv|| also contains the 
ABq^II site found for a metastable Mu~. A transition 
from Mu+ to Mu^ at these sites, observed above 150 
K, must involve a short-lived neutral that remains inside 
this cage. Related to the question of site, the origin for 
the high degree of anisotropy in the Mu hyperfine cou- 
pling is yet to be understood. A theoretical study of the 
shallow-donor N-H complex in diamond may prove 
helpful on this issue. 

More than 60% of implanted muons form diamagnetic 
states even at the lowest temperature. Our data indicate 
that whenever the shallow Mu*^ is present, the diamag- 
netic signal shows a fairly rapid exponential relaxation, 
characteristic of fast dynamics such as a transition out of 
that state. One possibility is that either Mu+ or Mu~, 
both of which were identified in the previous experiments 
[l3^|, represents the initial state formed upon implanta- 
tion and is a precursor to slightly delayed formation of 
the shallow neutral Mu center. 

Further experiments, including muonium level crossing 
resonance measurements, are clearly required to clarify 
the detailed electronic structure and to more precisely 
identify the site of the observed shallow Mu*^ center. 

In summary, we have demonstrated that a paramag- 



netic muonium center with extremely small hyperfine pa- 
rameters is formed in GaN below 25 K. This Mu° cen- 
ter has a hyperfine interaction that is axially symmet- 
ric around [0001]. The temperature dependence of the 
yield for this state indicates that it acts as a shallow 
donor, strongly suggesting that hydrogen, if located at 
an equivalent crystallographic site, would contribute to 
the unintentional n-type conductivity in GaN. 

We would like to thank the staff of TRIUMF for their 
technical support. We are also grateful to K. Nagamine 
and K. Nishiyama for their continuous encouragement 
and helpful discussions. This work was supported by the 
Natural Sciences and Engineering Researech Council of 
Canada (KHC, BH), the US National Science Foundation 
(RLL), and the Robert A. Welch Foundation (RLL). 



* Also at School of Mathematical and Physical Science, 
The Graduate University for Advanced Studies 

t Present address: Photonic and Wireless Devices Re- 
search Laboratories, NEC Corpration, Ohtsu, Shiga 520- 
0833, Japan 

* Present address: Computational Materials Science Cen- 
ter, National Institute for Materials Science, Tsukuba, 
Ibaraki 305-0047, Japan 

[1] H. Amano et al, Jpn. J. Appl. Phys. 28, L2112 (1989). 

[2] S. Nakamura et al., Jpn. J. Appl. Phys. 31, 1258 (1992). 

[3] For a recent review, see, for example, Ill-Nitride Semi- 
conductors: Electrical, Structual and Defects Properties, 
ed. by O. Manasreh (Elsevier, Amsterdam, 2000). 

[4] H.P. Maruska et al., Appl. Phys. Lett. 15, 327 (1969). 

[5] M. Ilegems et al., J. Phys. Chem. Solids 34, 885 (1973). 

[6] J. Neugebauer et al, Phys. Rev. B 50, 8067 (1994). 

[7] D.C. Look et al., Phys. Rev. Lett. 79, 2273 (1997). 

[8] W.J. Moore et al, Phys. Rev. B 65, 081201(R) (2002). 

[9] J.A. van Vechten et al, Jpn. J. Appl. Phys 31, L139 
(1992). 

[10] K.H. Chow, B. Hitti, R.F. Kiefl, in: Identification of 
Defects in Semiconductor, M. Stavola (Ed.), (Academic 
Press, New York, 1998). p. 137. 

[11] R.L. Lichti et al., Physica B 289-290, 542 (2000). 

[12] R.L. Lichti et al, Physica B 308-310, 73 (2001). 

[13] R.L. Lichti, Physica B 326, 139 (2003). 

[14] J.M. Gil et al., Phys. Rev. Lett. 83, 5294 (1999). 

[15] S.F.J. Cox et al., Phys. Rev. Lett. 86, 2601 (2001). 

[16] J.M. Gil et al., Phys. Rev. B 64, 075205 (2001). 

[17] K. Shimomura et al, Phys. Rev. Lett. 89, 255505 (2002). 

[18] A.F. Wright, Phys. Rev. B 60, R5101 (1999). 

[19] S.M. Myers et al, J. Appl. Phys. 88, 4676 (2000). 

[20] J. Neugebauer et al., Phys. Rev. Lett. 75, 4452 (1995). 

[21] C. G. Van de Walle et al.. Nature 423, 626 (2003). 

[22] T. Miyazaki et al, Phys. Rev. Lett. 88, 066402 (2002). 



